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Abstract 

We calculate exact analytical expressions for 0(a s ) 3-jet and 0(a 2 ) 
4-jet cross sections in polarized deep inelastic lepton nucleon scatter- 
ing. Introducing an invariant jet definition scheme, we present differ- 
ential distributions of 3- and 4-jet cross sections in the basic kine- 
matical variables x and W 2 as well as total jet cross sections and 
show their dependence on the chosen spin-dependent (polarized) par- 
ton distributions. Noticebly differences in the predictions are found for 
the two extreme choices, i.e. a large negative sea-quark density or a 
large positive gluon density. Therefore, it may be possible to discrimi- 
nate between different parametrizations of polarized parton densities, 
and hence between the different physical pictures of the proton spin 
underlying these parametrizations. 



I. INTRODUCTION 



The spin structure of the longitudinally polarized proton has attracted much 
theoretical interest during the past years. This is mainly due to the disagree- 
ment between the theoretical prediction of the Ellis- Jaffe sum rule [|]] for the first 
moment of the polarized structure function gi(x,Q 2 ) and the combined SLAC- 
EMC-data 0] . The experimental results which turn out to be much smaller than 
expected suggest two extreme scenarios: a large negative polarization for the 
sea-quarks [||] or alternatively a large positive polarization of gluons in the longi- 
tudinally polarized proton [|],[|] (see also the discussion in sect. II). A combination 
of both is also possible. 

There have been many suggestions for measuring polarized parton distribu- 
tions with hard processes . Only processes where the polarized gluon distribu- 
tion Ag enters at leading order are a promising source for measuring the polarized 
parton distributions. 

In this paper we study 3- and 4-jet cross sections in polarized deep inelas- 
tic lepton-nucleon processes using recently proposed polarized quark and gluon 
distributions, Ag and Ag respectively. 

To 0{a s ) one has 3-jet final statesB. The relevant subprocesses are 7*g — > 
Gq and 7*G — > qq. In 0(a 2 s ) the following subprocesses contribute to the 4- 
jet production: 7*g — ► GGq, r y*q — > qqq and 7*G — ► qqG. All partons are 
assumed to be massless. We present exact analytical expressions for the relevant 
3- and 4-jet partonic hard scattering cross sections, using polarized beams and 
polarized targets, and specify how these partonic cross sections are folded with 
the respective polarized parton densities to obtain the DIS hadronic cross section. 
Introducing an invariant jet definition scheme we present differential distributions 
of 3- and 4-jet rates. 

Large differences in the differential distributions are found for the two ex- 
treme choises for the polarized parton distributions made in the literature i.e. a 
large negative sea-quark density or a large positive gluon density. Thus, the anal- 
ysis of jet production in polarized deep inelastic scattering allows to extract Aq 
and/or AG from future measurements and it may be possible to discriminate 
between the distinctly different physical pictures of the proton spin underlying 
these parametrizations. 

The paper is organized as follows: In section II we introduce our notations 
and some useful formulae. In section III we discuss the problem of defining jets in 
deep inelastic scattering and present analytical results for the 0(a s ) 3-jet cross 



*We include the target jet when counting the number of produced jets. 
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sections. Section IV contains analytical formulae for the 0(a 2 ) 4- jet cross sections. 
We specify in detail how the partonic hard scattering cross sections are folded 
with the respective parton distribution functions. In setion V the set of polarized 
parton distribution functions used in our calculations are reviewed. Section VI 
contains our numerical results. We discuss overall kinematical cuts that we choose 
to define our deep inelastic scattering phase space. Further cuts are then imposed 
on the invariant jet-jet masses to define the resolved 3- and 4-jet phase space. 
We present predictions for the 3- and 4-jet cross sections and show distributions in 
the kinematical variables x and W 2 . Section VII contains our summary. Finally, 
detailed analytical formulae for the 4-jet partonic matrix elements are relegated 
to the appendix. 



II. KINEMATICS OF POLARIZED DIS AND THE STRUCTURE 

FUNCTION Gi 

Consider deep inelastic lepton-proton scattering (X is an arbitrary hadronic 
final state) 

l 1 (l)+H(P)^l 2 (l')+X(P f ) (1) 

The particle momenta are given within the brackets. 

Reaction (|I|) proceeds via the exchange of an intermediate vector boson V = 
7*, Z, W^. In this paper we will only consider the exchange by a virtual photon. 
We denote the 7*-momentum by q, the absolute square by Q 2 , the cms energy by 
s, the square of the final hadronic mass by W 2 and introduce the scaling variables 
x and y: 



q = 


l-V 


Q 2 = 


—q 2 = xys > 


s = 


{P + lf 


W 2 = 


Pf = (P + q) 2 


X = 


o 2 


y = 


§ (0 < V < i). 



For fixed s, only two variables in (0) are independent, since e. g. 

xW 2 = (1 - x)Q 2 , Q 2 = xys. 

For polarized lepton and polarized target scattering, the inclusive (with respect 
to the outgoing hadrons) cross section for process (0) can be written as 
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dAa dV T d 2 a n 



,2 



■Sna" f Mxy] , ^ 9n 2Mx , 

\2-y- -^-\gi{x, Q) - —92(x, Q l 



dx dy dx dy dx dy q 2 



(3) 



The left arrow above a denotes the polarization of the incoming lepton with 
respect to the direction of its momentum. The right arrow stands for the polar- 
ization of the proton parallel or anti-parallel to the polarization of the incoming 
lepton. gi(x,Q 2 ) and g2(x,Q 2 ) are the polarized structure functions for longitu- 
dinal and transverse spin respectively. At large energies, only gi contributes. 
In the QCD improved parton model, the polarized structure function gf of the 
proton can be expressed in the lowest-order approach byi 

2g{ (x, Q 2 ) = ^Au v (x, Q 2 ) + ^Ad v (x, Q 2 ) + ^Aq(x, Q 2 ) (4) 

We have taken an SU(3)-flavour symmetric polarized sea Aq = Au = Ad = As = 
As where Au v and Ad v are the polarized valence quark distributions and 

Aq(x,Q 2 )=q + (x,Q 2 )-q4x,Q 2 ). (5) 

In (0) q+(q~) denotes a parton of flavour q in the proton whose spin is aligned 
(anti-aligned) to the proton's spin. 
Experimental measurements of the asymmetry 

provide us with a determination of gf (a;, Q 2 ). In (|]) F^(x, Q 2 ) is the unpolarized 
structure function, i.e. F<[(x, Q 2 ) = xJ2 e l[q( x , Q 2 ) + <l( x i Q 2 )}- 

1 

Writing gf (Q 2 ) = J g[(x,Q 2 )dx as 
o 

9i(Q 2 ) = \\\^{Q 2 ) + l -Ad{Q 2 ) + l -As{Q 2 )] (7) 

and using the F/D ratio of the hyperon decay and the Bjorken sum rule, the 
following values for the first moments of the polarized quark densities have been 
obtained 

Am = 0.78 ±0.06 
Ad = -0.47 ±0.06 
As = -0.19 ±0.06 



^We follow the notation of Gliick, Reya, Vogelsang 
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and 



Ag = Au + Ad + As = 0.12 ±0.17. (8) 

Therefore the total spin carried by the quarks and antiquarks in a polarized 
proton is small and actually compatible with zero. This is the famous EMC result. 



Recently, new data from the SMC muon-deuteron scattering experiment pi and 



from the SLAC muon-He scattering experiment |l Ij are available. It is shown in 
Yl\, that a global average of the present existing date yields Ag = Au+Ad+As = 
0.22 ± 0.10 which is still not very different from the value in (|5|). 

Several groups have suggested that the difference between the Ellis- Jaffe sum 
rule |l| (the sea contribution was set to zero in the Ellis Jaffe sum rule) and 
the EMC result could be resolved by a larger higher correction due to a large 
gluon polarization or by a large sea contribution ||. Various sets of helicity 
difference distribution functions have been developed based on these ideas. An 
independent determination of the gluon spin densities are therefore crucial in the 
understanding of the spin structure of the proton. 

In the following, we study 0(a s ) 3-jet and 0(a 2 s ) 4-jet production in polarized 



DIS (using the polarized parton distributions proposed in JT3[) and show that 
these processes yield strong information on the spin-dependent quark and gluon 
densities. 



III. JET CROSS SECTIONS IN 0{a s ) 
The Quark-Parton model (O (a°)) process 

h(l) + q(p = r)P) — fe(0+?(Pi) 

predicts two and only two jets in the final state, namely the struck quark jet and 
the hadron remnant jet which leads to the prediction of two back-to-back jets in 
the hadronic center-of- mass-frame. 

The 0(a s ) tree graph processes (see fig. la,b for the relevant boson-parton 
subprocesses) 

h(l)+ V (P) — k(l')+ V (ft) + G(p 2 ) (9) 
h(l) + G(p) — > l 2 (l') + q( Pl ) + q(p 2 ) (10) 

give rise to 3-jet production. There are singular regions in the phase space of 
the 0(a s ) tree graph cross sections. The singular behaviour corresponds to the 
situation where two (or more) partons are irresolvable and have to be counted as 
contributions to the 2-jet cross section. Parts of these singular cross section con- 
tributions will cancel against corresponding singular terms of the 0(a s ) virtual 
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2-parton corrections. A remaining initial state mass singularity has to be absorbed 
into the polarized parton densities. 

In order to calculate a 3-jet cross section we have to define what we call 3 
jets by introducing a resolution criterion. As has been elaborated in detail in 
IP1,T5[] energy- angle cuts are not suitable for an asymmetric machine with its 



strong boosts from the hadronic cms to the laboratory frame. As a jet resolution 
criterion we use the invariant mass cut criterion introduced in [1T4| such that 



% > M 2 C = max{y cut W 2 , M 2 } = max{y cut { Pl + p 2 + p 3 ) 2 , M 2 } = 1, 2, 3;i # j) 

(11) 

where y cut is the resolution parameter and is the invariant mass of any two 
final state partons, including the remnant jet with momentum P3 = (1 — rj)P. 
Mq is a fixed mass cut which we have introduced in order to clearly separate 
the perturbative and non-perturbative regime in the case where W 2 is small (see 
the detailed discussion in sect. VI). The 3-jet phase space is represented by the 
kinematical variables 

Q 2 x ppi 

x„ = = — z — and <p. (12) 

2 pq 7] pq 

The angle denotes the angle between the parton plane (p,pi) and the lepton 
plane (1%, l 2 ) in the p + q = center of mass frame. 

The 0(a s ) 3-jet cross section for polarized DIS is then obtained by 

d 2 Aa _ 27ia 2 (2 - y) a s {fi 2 R ) V J X ' *7" dz J <ty ^ 



dxdy Q 2 2n J r] J J 2ir 

Vmin z min 



{/S + ^+^cos^(/P» + /«)}. 



r (o) , r (o) , Vi -_y 
y 

The integration limits in 77 and z originate from our jet resolution criterion ([II 
such that 

M? 



■~iinii , -, \ 1 Zmaxi ( 14 ) 

ys[l -77) 

Ml 2M C 2 
ys ys 

The 3-jet integrands are given by (i = 0, 1) 

n f 

J S = E e /( A ^,^) + Aq( v ,fi 2 F ))A\M^ qG , (16) 
/=i 

n f 

^ = (E4) AG (^^) A|-M«| 2 G ^ (17) 
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where fip denotes the factorization mass in the polarized parton densities Aq and 
AG. The 3-jet matrix elements A\Ai^\g^ qg and A\Ai^\Q^ q g are given by 



x 2 + z 2 



A|^ (0) |^ G = C F [ {1 _ P Xp){1 _ z) + 2(x p + z)}, (It 



A|M (0) |^ W - = T R (2x p - l)(- z + ^-2), (19) 



*\M^ = 4CJ £ (1-,,-z), (20) 



A\M^\ 2 G ^ = 4T R 



x p (l x p ) 



(21) 



where Cf = I and Tr = | are the corresponding colour factors. Note that the 
contribution of the matrix elements AIA^ 1 - 1 ] 2 vanish after integration of the angle 
in (|T3|) . Furthermore, the z- integration can be done analytically yielding: 



/n 27T ig^g^ * i _ 1 _ 

3 1 

2 1 - x 



+ \2x p + l-\—^— ) (1-2^)1, (22) 



1 ^A|M(°)|^ W -^ = -T fi (2x p - 1) (in + 1 - 2z miT . 



IV. FOUR JET CROSS SECTIONS 

In this section we provide exact 0(a 2 ) 4-jet production cross sections in polar- 
ized deep inelastic scattering. The lowest order contributions to the 2 — > 4 parton 
process Z + p — >■ /' + pi + ^2 + £>3 are (the relevant boson-parton subprocesses are 
depicted in figs. (2a-c)): 

l + q^l' + q + G + G, (23) 
l + G—>l' + q + q + G, (24) 
l + q — >l' + q + q + q (25) 

and antiquark initiated processes. The corresponding differential partonic cross 
sections dAa p are given by (s p = (p + I) 2 ) 



= i A^p! dPS «, 
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Here A\A4^~-' et \ 2 denotes the squared matrix element (without the photonprop- 
agator) of the parton processes (|23"D - (p5|) including all coupling, statistical and 
colour factors. The lorentz invariant four particle phase space 

dPsi*> = (2K)W(!+p-v- Pl -p 3 -p a y d3f ' d ^ 1 d ^ 2 d ^ 



(2ir) 3 2E' (2tt) 3 2E 1 (2tt) 3 2£ 2 (2it) 3 2E 3 



can be parametrized by seven kinematical variables |15|], which we choose as 
y, Q 2 , z', (ft, s 23 , cos 6*23 and 2 3 where y and Q 2 are already defined in <^). s 23 
denotes the squared invariant mass of partons 2 and 3 

S23 = (P2 +P-aY 

and the variable z' is given by 

/ = 1 PPi 
pq ' 

4> is the azimuthal angle of the outgoing lepton in the (boson-initial parton) center 
of mass frame. Finally 8 23 and 023 denote the polar and azimuthal angles of parton 
2 in the rest frame p 2 + f>3 = and s — (p + q) 2 . In terms of these variables the 
phase space can be decomposed as [T3|] 

1 1 ys 1 1 27T z's +1 2tt 

jd-qj dPS& = J dy JdQ 2 Jd V J dz' Jd<j) J ds 2 3 Jd COS #23 / #23- (26) 

OOzOOO-1 

The 0(a 2 ) hadronic 4-jet cross section is then given by 



(27r) 2 d*Aa H g , q.( A «g t ) ^7ra a 2 A[^ e ^ 

dydQ 2 dr/dz'd(f)ds 2 3d cos 6 2 3d<p 2 3 i ^f n ^ 2tt ' 8r/s % ' Q 4 

(27) 

where the index i refers to the nature of the incoming polarized parton (gluon: 
i — 0, rif light quarks: i = 1, . . . rif and rif light antiquarks: i = —rif, . . . — 1). 
fip denotes the factorization mass in the polarized parton densities and [xr is the 
renormalization mass entering the strong coupling constant. 
In the next section we explicitly specify the squared matrix elements A| i M^ _: ' e *| 2 
including their respective parton weights fi(r), jji 2 F ) for the partonic processes (p3|)- 



A. 0(a 2 ) Four jet matrix elements in polarized eP scattering 

We are now in the position to calculate the squared matrix elements of the 
hard scattering process. In order to exhibit the parton that initiates the hard 
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scattering process we write down the relevant hard scattering process with the 
parton weights included. We will also explicitly include the colour factors and 
final state identical particle factors in our expressions. In the following we specify 
the various quark, antiquark and gluon initiated processes: 



(i) Quark and antiquark initiated two quark two gluon processes 

(<?/ — * qfGG) and (q~f — > q/GG) 

These processes receive contributions from the diagrams of fig. 2a. We divide the 
contributions into four classes, \-M-\ ABClC2 distinguished by their group weight 
factors (we follow the notation in f|T7| for e + e~ jet production): 

A : planar abelian graphs: NcCp 11, 21, 22, 31, 32, 33, 

44, 54, 55, 64, 65, 66 

B : non-planar abelian graphs: N C C F (C F - ) 41, 42, 43, 51, 52, 53 

61, 62, 63 

G x : QED-QCD interference graphs: \N C C F 71, 72, 73, 74, 75, 76, 

81, 82, 83, 84, 85, 86 
C 2 : QCD-QCD graphs: N C C F 77, 87, 88. 

The squared matrix element is then given by: 
1 1 



^\M%t qfGO ? = ;T q 5>/ + A ?/) I ® A \M qf -, qf GG\ 2 A+B+Cl+ c 2 (28) 



The numerical factors in (28) are initial state averaging factor for parton colour 



and the final state identical particle factor. The symbol (g> denotes the folding 
with the parton densities. The numbers ij in the third column refer to the inter- 
ference of graph i with j (= > j) in fig. 2a. The explicit expressions for the 
squared matrix elements ij can be found in the Appendix. They are denoted by 
d(i,j) and include the color, statistical and average factors. 

(ii) Gluon initiated process (G — > Gqfq~f) 

These processes receive contributions from the diagrams of fig. 2b. The contribu- 
tions can be divided in similar classes as before (Colour factors and associated 
graphs are explicitely given): 

A : planar abelian graphs: 11, 21, 22, 31, 32, 33, 

TrC f (Nq - 1) 44, 54, 55, 64, 65, 66 

B : non-planar abelian graphs: 41, 42, 43, 51, 52, 53 

T R {C F -^){Nl-l) 61,62,63 

C x : QED-QCD interference graphs: 71, 72, 73, 74, 75, 76, 

\T R N c (Nl - 1) 81, 82, 83, 84, 85, 86 

C 2 : QCD-QCD graphs: T R N C (N% - 1) 77, 87, 88. 
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The squared matrix element reads 



A|A<g% ff/ 1 2 = g • I E 4 I AG ® ^^lW 1+ c 2 (29) 




where | is an averaging factor for the initial gluon colours. The numbers ij in 
the third column refer to the interference of graph % with graph j in fig. 2b. The 
explicit expressions for the squared matrix elements ij are to long to be presented 
in this paper. They are available from the authors. 

(iii) Quark and antiquark initiated four-quark process (q f — > qfqpqp) and 
(q f -> q f qpqp) 

These diagrams receive contributions from the diagrams in fig. 2c. The contribu- 
tions are divided in four classes distinguished by their group weight factors: 

D : non-interference contributions: N C C F T R 11, 21, 22, 55, 65, 66 
D' : non-interference contributions: N C C F T R 33, 43, 44, 77, 87, 88 
E : non-interference contributions: 31, 32, 41, 42, 51, 52, 

N C C F (C F - ^p) 61, 62, 73, 74, 75, 76, 

83, 84, 85, 86 

F : interference contributions: NcC f Tr 53, 54, 63, 64, 71, 72, 



81, 82. 



The squared matrix element reads 

1 lr»/ 



' I q^qqq 3 2 



J2e}(Aq f + Aq f )®(np.\M\l + \M\l) 



+ E(E4)( a ^/ + a ^/)®I- m Id' (so) 
/=i /'=i 

+ E ( E e /')e/ (A 9/ - A^) ® (Ml 2 , 
/=i /'=i J 

(| colour averaging, | statistical factor for identical particles in the final state). 
The numbers ij in the third column refer to the interference of graph % with 
graph j in fig. 2c. The explicit expressions for the squared matrix elements ij 
can be found in the Appendix. They are denoted by e(i,j) and include the color, 
statistical and average factors. 

Some comments are in order for this subprocess: We do not distinguish be- 
tween equal and unequal flavour cases. The overcounting in the unequal flavour 
case (/ 7^ /') is compensated by our multiplication of these contribution with the 
Fermi statistical factor \ as in the equal flavour case (/ = /'). 
The interference contributions class F involve the calculation of two fermion 
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traces. They involve a product of traces of 3 fermion (or antifermion) propaga- 
tors. Each such trace is anti-symmetric under quark-antiquark exchange for the 
vector current 7^ due to charge conjugation invariance. Note that the interference 
class F does not contribute to the total cross section if one does not identify the 
flavours of the produced quarks and antiquarks. This is a consequence of charge 
conjugation - the non-abelian generalization of Furry 's theorem. 



V. POLARIZED PARTON DENSITIES 

The hadronic 3- and 4-jet cross sections are obtained by folding the partonic 



cross sections with the polarized parton densities (see ( |T3| . f2~7D ) . In this section, we 
review the set of the polarized parton density parametrization |13| which we will 
use in our numerical calculation. The sea and gluon parton densities are chosen in 
such a way that the experimental result for the Ellis- Jaffe sum rule is reproduced 
up to 0(a s ) . 

These densities are given by: 

Au v (x) = a(x)u v (x,M 2 ) (31) 
Ad v (x) = {3(x)d v (x,M 2 ) (32) 



with 



a(x) = x - 287 , (3(x) = (?-^)z 

v 1 - xo 



p 



with xq = 0.75 and p = 0.76. The M 2 ) and dy(x, M 2 ) are the unpolarized 

parton densities given by the DFLM4 set (A=0.2 GeV) at M 2 = lOGeV 2 |TE 
According to |13|] we can adopt the following parametrization: 



set 1: At M 2 = 10 GeV 2 one has a large negative polarized sea-quark density and 
a zero polarized gluon density 

As(x) = -11.8 x a94 (l - x) 5 s(x, M 2 ), 
AG(x) = 0. (33) 

set 2: At M 2 = 10 GeV 2 one has a zero polarized sea-quark density and a large 
positive gluon density 

As(x) = 0, 

AG(x) = 6.0 x a76 (l - x) 3 G(x, M 2 ). (34) 
The above gluon density is about the same as the one proposed in |5[]. 
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set 3: In a realistic situation it is unlikely that As(x) or AG(x) vanishes at some 
scale for all x. Therefore the authors in |l3f also proposed 

As(x) = -3.39 a; - 62 (l - x) 1A s(x, M 2 ), 
AG(x) = 2.69 x°- 76 (l - x) 3 G(x, M 2 ), (35) 

for M 2 = 10 GeV 2 . 

s(x, M 2 ) and G(x, M 2 ) represent the unpolarized parton densities belonging to 
the DFLM4 set @. 

VI. NUMERICAL RESULTS 

We will now turn to our numerical cross section results. In order to define the 
deep inelastic scattering region, we choose the lower bounds 

Q 2 > Ql = 4 GeV 2 , 

W 2 > Wl = 5 GeV 2 , (36) 
x > x = 10~ 2 . 

We introduce a limit in W 2 to insure an appropriate hadronic final state. When 
producing well separated multijet final states the effective lower limit of W 2 
will actually be increased (see ( |3*9| , fi0| )). As already discussed in section III one 
encounters infrared/mass (IR/M) singularities when calculating the 0(a s ) and 
0(a 2 ) tree graph cross sections, which originates from the emission of soft and/or 
collinear partons. The tree graph cross sections make sense only for the production 
of well separated hard jets. We introduce the following invariant mass-cut-off to 
define the finite tree graph cross sections! (see fllTD): 

Sij > max{y c W 2 , M 2 }, Sy = (p< + p 3 ) 2 i^j = 1,2, 3, (4) (37) 

where y c is the resolution parameter for which we choose y c = 0.04 [pj4] , |T5| . Mo 
is an additional fixed mass cut which we have introduced in order to clearly 
separate the perturbative and non-perturbative regime in the case, where W 2 
is small. Note that W 2 is limited from below by Wq (see fl3~6"D ). This leads to 
too small invariant masses for a perturbative calculation to be meaningful. For 
example with y c = 0.04 one obtains < 0.2 GeV 2 . A reasonable choice for the 



*Note, that the remnant jet (p 3 in 0(a s ) and p 4 in 0(a 2 ) ) is included in this 
definition. 
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fixed mass cut in fl3~T|) is M = 2 GeV [15[]. Our subsequent analysis is then based 
on the invariant mass cut 



Sij > max{0.04 W 2 , 4 GeV 2 }. (38) 

The invariant mass cut in (^8|) introduces then the following available ranges in 
rj and W 2 for the 0(a s ) 3-jet and 0(a 2 ) 4-jet cross sections: 
0(a s ) 3-jet case: 

W 2 > max{Wl 3 M 2 }= 12 GeV 2 
V >y c + x(l-y c ) (39) 



0(a 2 s ) 4-jet case: 



W 2 > m&x{W$, 6 M 2 } = 24 GeV 2 
V> Sy c + x(l-3y c ) (40) 



As can be seen from ( |39"D and (0), the lower bound in rj and W increases with 
jet multiplicity. 

In our numerical calculation we use the one-loop formula of the strong coupling 
constant 

2 127T 

= : — 7x 



(33 -2n f ) ln^ 

and the A value consistent to our choise of polarized parton distributions. The 
renormalization scale fi 2 R and the factorization scale fi 2 ? are both fixed to Q 2 . 
The value of a s is matched at the thresholds q = m q and the number of flavours 
rif in a s is fixed by the number of flavours for which the masses are less than 
Q 2 . Furthermore the number of quark flavours that can be pair-produced are set 
equal to rif chosen in a s . 

In figs. 3 a-f we show how the various cross sections vary with a given invariant 
energy yfs for the three different sets of parametrizations (set 1-3) of the polarized 
parton densities. The Q 2 range is chosen as 5 GeV 2 < Q 2 < s. One observes, 
that all quark initiated subprocesses are positiv and have a maximum around 
yfs = 15 — 20 GeV. In contrast, the gluon initiated contribution to Act is negativ 
and decreasing with increasing y/s. 

In figs. 4-7 we show the dependence of the 0(a s ) 3-jet cross sections on the 
basic kinematical variables x and W 2 for the different parametrizations. We have 
restricted the Q 2 range from 5 GeV 2 < Q 2 < 7 GeV 2 (figs. 4,6) and 9 GeV 2 < 
Q 2 < 11 GeV 2 (figs. 5,7). The quark (dashed) and gluon (dotted) initiated cross 
sections are shown separately as well as the sum (solid). 



12 



Comparing figs. 4a-c (5a-c) one observes large differences in the x distribution 
using different parametrizations. In particular, figs. 4b (5b) show that the 3-jet 
cross section (solid curve) is negative for a large gluon contribution for x values 
less than about 0.04. An analysis of the x-dependence at lower x values is therefore 
very sensitive to the polarized gluon densities. 

Figs. 6a-c (7a-c) show the lU-distribution of 3-jet production for the three 
different sets of parametrizations, where the Q 2 range is given by 5 GeV 2 < Q 2 < 
7 GeV 2 (9 GeV 2 < Q 2 < 11 GeV 2 ). Here larger values of W are more sensitive to 
the choice of polarized parton distributions. The gluon process gives a negative 
contribution and dominates the resulting cross section in the large W range in 
figs 6b, c. This is a reflection of the x behaviour since W 2 = (l—x)ys. Like for the 
^-distribution, the predictions for the two extreme choices of parton distributions 
are very different. 

Figs. 8 and 9 show predictions of the x and W dependence of the 0(a 2 ) 4- 
jet cross sections for the three different sets of parametrizations and 5 GeV 2 < 
Q 2 < 7 GeV 2 . The behaviour of the x and W dependence is similar to the 0{a s ) 
3-jet cross section. One observes that the q — > qGG subprocess dominates the 4 
quark process q — > qqq over the whole x and W ranges. However, for set 2 of the 
parton distributions (fig. 8b), the gluon process G — > qqG gives a large negative 
contribution for smaller x values, which gives again a negativ cross section for 
the sum of all contributions. The sum of the contributions for set 3 of the parton 
distributions (figs. 8c, 9c) remains still positive indicating that the 4-jet cross 
sections are less sensitive to the different parametrizations. The behaviour of the 
W dependence in the 4-jet distributions is again similar to the behaviour found 
for the 3-jet cross setcions. The gluon process gives a negative contribution at 
larger W values. 

VII. CONCLUSIONS 

Jet production in DIS scattering experiments, using polarized beams and po- 
larized targets, has been studied. It is shown that the jet cross sections are strongly 
dependent on the particular form of the spin dependent parton distributions used 
in the calculation. Therefore, jet production in polarized DIS can be used to put 
rigid constraints on the various parametrizations of polarized structure functions. 
These informations allow to discriminate between the different physical pictures 
of the proton spin underlying these parametrizations. 
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APPENDIX A: MATRIX ELEMENTS 



We give here explicit formulae for the squared matrix elements A\A4^ ■ ?et | 2 in 
(p^j- p0|). As usual, A| < M i ~"' e *| 2 can be written in terms of the lepton tensor L^ LU 
and the hadron tensor AH^ = \{H^ V + - ~), where + (H> XV ~) denotes 
the hadron tensor for the process when the incoming parton has positiv (negativ) 
helicity. The hadron tensor AH^ V can be expanded in terms of the following 
tensors: 



[p,Pi] 

\P,P2] 

\p,Pa] 

[Pl,P2] 

IPhPs] 

\P2:Ps] 



e^ al3 p a p 2 /3 



(Al) 



Furthermore, we introduce the following notation 



s 


■= PP3 


t 


■= P2P3 


u 


:= pp 2 


X 


■= P1P3 


y 


■= pp 1 


z 


■= P1P2 



(A2) 



To obtain the full matrix elements A\A4 



A- jet 1 2 



the corresponding hadron tensors 



for each subprocess have to be contracted with the anti-symmetric part of the 
the lepton tensor L^ v : 



[IV 



2ie 



fiVOlfl 



(A3) 



1. Matrix elements for (<?/—» qfGG) and (<?/—» qfGG) 

The results for the interference of graph i with graph j in fig. 2a is denoted 
by d(i,j). Note that d(i,j) contain the color, statistical and average factors as 
given in (^8|) . 

d(l,l) = -16ifc> 2 ,p 3 ]/(9sz) 

d(2, 1) = 16i(\p,px]u+ \p,p 2 ](s - x + y + 2u - t) + [p,p 3 }u 

+ [pi,p 3 }u + [p 2 ,p 3 ]uj/(9sz(x + z + t)) 
d(2, 2) = -16i( - \p, Pl ]{z + t) + [p,p 2 }x - \p,p 3 ](z + t))/(9z(x + z + tf) 
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d(3, 1) = 16i(\p,p 1 ]x+ [p,p 2 ]x+ \p 1 ,p 2 \x + \pi,p 3 ](2x - y + z - u + t) 
+ [P2, Pz]x) I (9sz(s + u-t)) 

d(3, 2) = lGi([p,pi](sy + su + x 2 — 3xy + xz — 2xu + 2xt + 2y 2 — yz + 3yu 
—2yt + u 2 - 2ut + t 2 ) + [p, p 2 ] (—2xy - xu + y 2 - yz + yu) 
+\p,P3](-xy + y 2 + yu - yt + zu - zt) + \p 1 ,p 2 \(sx + st - xy 
+y 2 + yu- yt) + [pi,p 3 ](sy - xy - xu + y 2 + 2yu) 
+ [P2,P3](sy + sz - xy - xu - yz + yufj/(9sz(s + u — t)(x + z + t)) 

d(3,3) = 16i([p,pi](s-*) + [pi,p 2 }(s-t) - [ Pl ,p 3 ]u)/(9s(s + u-t) 2 ) 

d(A, 1) = -2i([p,pi]y(s -x + 2y-z + u) + \p,p^y{s -x + y) 

+[p,Ps}y(y - z + u) + [pi,p 2 }y(s -x + y) 
+\pi,P3]y(y -z + u)^/(9sxzu) 

d(4, 2) = 2i(\p,pi](sz + zu + ut) + [p,p 2 ](sz + xz — xu + yz) 

+ [p,p 3 ](sz - 2xy + xz - xu - yz - yt) + \p 1 ,p 2 \(-xu + yt) 
-\p\,Ps]xu + [p 2 ,p 3 }y(x + z)^/(9xzu(x + z + t)) 

d(4,3) = 2i([p,p 1 }(sx + sz- xt) - [p,p 2 ]xu + [p,p 3 ](-xu + yt) 

+\pi,P2](sy + sz + su - xu + 2yu - yt) + \pi,p 3 ](-sy + sz 
+su — xu) — [p 2 ,p 3 ]y{s + ufj /(9sxu(s + u — t)) 

d(5,2) =4i(- [p,pi]t+ [p,p 2 ]{x + z) + [p,p 3 ]{x + z))/(9xz(x + z + t)) 
d(5, 3) = — 2i([p, pi](s 2 — 2sx + 3sy — sz + su — 2st + x 2 — 3xy + xz 
-xu + 2xt + 2y 2 -yz + yu- 2yt) + [p, p 2 ]y(s -x + y) 
+ [p,P3]y(s -x + y) + \p 1: p 2 ]y(s -x + y) 
+\pi,P3\y(s -x + yf)/(9sx(x + z + t)(s + u - t)) 

d{6,3) = 4i([p, Pl ]t + [ Pl ,p 2 ]{s + u) + [ Pl ,p 3 ]{s + u))/{9su(s + u-t)) 
d(7, 1) = i(\p,pi\(2sx + sz + st — x 2 + 2xy — 2xz + xu — xt + Ayz + 2yt 

+2zu — zt — ut) + [p, p 2 ] (sx + 2sz — x 2 + xy — xz + 2xu — xt 

+2yz + yt + 4zu — 2zt) + [p, p 3 ] (—sz + 2xy + 2xz 

+2xu + yz + 2zu) + 2[pi,p 2 ]s(x + t) 

p 3 ] (—2sz + xy — xu + 2yz + 3yt + zu — ut) 

+ [P2,Ps](yz + 2yt + 2zu))/(szt(x + z + t)) 

d(7, 2) = 2i([p,p 1 ](-xz -xt + z 2 - 2zt - t 2 ) + [p,p 2 ](x 2 + xt + z 2 - zt) 
+2[p,p 3 ](x 2 + 2xz + xt + 2z 2 + zt)) /(zt(x + z + i) 2 ) 



16 



d(7, 3) = — i([p,Pi](s 2 — sx + 2sy + 2sz + 3xu — 2yu — 2yt + 2zu — zt 
—u 2 + Aut — 2t 2 ) + [p,p 2 ](sy + 2xu — yu — 3yt) 
+ [p,p 3 ](sy — yu + yt — 2zu + 2zt) + [pi,p 2 ](—3sx — — 3st 
—2yu + 2yt) + [pi,P3](— sy + 3xu — 3yu — yt + zu + -ut) 
+[p2,P3](— — 2s,2 + 2xu — yu — 3ytfj / (st(s + u — t)(x + z + 1)) 

d(7, 7) = -i([p,pi](2a; 2 + Uxz + 9xt + Az 2 + Wzt - t 2 ) + [p,p 2 ](-5x 2 + 3xz 
+3xt + 8zt) + [p,p 3 ](5xz + 8xt - 3z 2 + bzt)} / (t 2 (x + z + if) 

d(8, 1) = — i([p,pi](2sx — 4sy + sz — 2su + st + xu + xt — 2yu + 2yt + 2zu 
—zt — u 2 + ui) + [p, p 2 ] (sx — 2sy — 2sz — xu + xt — yu + 3yt) 
+2\p,p 3 ]z(u — t) + [pi,p 2 ](2sx — sy — sz + 2sw + 2xu — 2yu) 
+ [Pi ) Pz\ (^ sx ~ 2sy + 2sz — su + 2st + 2xu — yu + yt + zu — u 2 + ui) 
+ [p 2 ,p 3 }(2sx -sy + 2yt))/(szt(s + u - t)) 

d(8, 2) = i([p,pi](2sx + 2sz + st-x 2 + 2zy + 3ra - 4xt - 2yz - 2yt + z 2 
—zu - 2t 2 ) + [p, p 2 ] (sx + 3xy + 3xu - xt + yz - yt) + [p, p 3 ] (-sz 
+2xy + 2yt — 3zu + 3zt) + [pi,p 2 ](— 2sx — 2st + xy — yz + yt) 
+ [Pi,P3}(xy + 2xu - yz - 3yt) + [p 2 ,p 3 ](-2sz + xy + 2xu + yz - 3ytfj 
■1 / (zt(s + u - t)(x + z + t)) 

d(8, 3) = -2i([p,p!](s 2 - su + 2st + ut - t 2 ) + 2[p u p 2 ](2s 2 + 2su - st + u 2 - ut) 

+ [Pi,P3](s 2 + st + u 2 - ut))/(st(s + u- t) 2 ) 
d(8, 7) = -i([p,pi](Asx + Usz + st + lAxu - xt - lQyt + zt - ut + 2t 2 ) 

+ \p,p 2 ](—5sx + 7xu + 5xt — 8yt) + [p,p 3 ](5sz — 8yt — 7zu + 7zt) 

+ \pi,p 2 ](—5sx + 7sz — 5st — 8yt) + [pi, p 3 ](5xu — 8yt — 7zu — 7ut) 

+ 12[p 2 ,p 3 }(-sz + xu))/(t 2 (s + u- t)(x + z + t)) 
d(8,8) = -i([p,pi](2s 2 + Usu-9st + 4u 2 - llut-t 2 ) 

+ [pi,p 2 ](—5s 2 + 3su — 3st — 8ut) + [pi,P3](5su — 8st — 3u 2 — 5ut)^j 

■l/(t 2 (s + u-t) 2 ). 

The remaining matrix elements can be obtained by the interchange of 
P2 P3, s <-> u, x <-> z from the above results. 



d(4,4) 


- d(l,l) 


d(5,4) <- 


— > d(2,l) 


d(5,5) <- 


-»• d(2,2) 


d(6,4) <- 


-»• d(3,l) 


d(6,5) <- 


— »• d(3,2) 


d(6,6) <- 


-> d(3,3) 


d(5,l) <- 


- d(4,2) 


d(6,l) <- 


— > d(4,3) 


d(6,2) <- 


-> d(5,3) 


d(7,4) - 


- d(7, 1) 


d(7,5) - 


d(7,2) 


d(7,6) <- 


-> d(7,3) 


d(8,4) <- 


- d(8,l) 


d(8,5) <- 


— > d(8,2) 


d(8,6) <- 


-> d(8,3) 



(A4) 
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2. Matrix elements for (G — > Gqjqj) 

The result for the interference of graph z with graph j in fig. 2b is denoted by 
c(i,j). The analytical results are to long to present here. They are available from 
the authors. 

3. Matrix elements for (qj — > QfQfQf) and (qy — ► qtqfig.fi) 

The results for the interference of graph i with graph j in fig. 2c is denoted by 
e(i,j). Note that e(i,j) contain the color, statistical and average factors as given 
in ©. 

e(l, 1) = — 2i(\p, Pi](— 2xz — xt — zt) + [p, p 2 ](x 2 — xz — zt) 

+ [p,p 3 ](-xz -xt + z 2 ))/(3t 2 (x + z + t) 2 ) 
e(2, 1) = 2i(2\p, Pi](— sz — xu + yt) + [p,p2}(sx — xt — xu + yt) 

+\p,P3](— sz +yt — zt + zu) + \p 1 ,p 2 ](sx — sz + st + yt) 

+ [Pi,P3](-xu + yt + zu + tu) 

+2\p 2 ,p 3 ](sz - xu))/(3t 2 (x + z + t)(s-t + u)) 

e(2, 2) = — 2i(\p,pi](st — 2su + tu) + [pi,P2\(s 2 — su + tu) 

+ [Pi,P3](st -su + u 2 ))/(3t 2 (s - t + uf) 
e(3, 1) = 2i(\p,pi]z(x + t) + [p,p 2 ]z(x + t) - \pi,p 3 ]z(y + u) 

-\P2,P3\z(y + u)J/(9tu(x + z + t)(y - z + u)) 
e(3,2) = 2i(\p,p 1 ](-xu + yt) + \p,p 2 ](-sz - xu + yt) + [p,p 3 }z(-t + u) 

+\pi,P2]s(y + u) + [pi,p 3 }(sz -xu + yt) 

+ [P2,P3]{~xu + yt)^/(9tu(s -t + u)(y -z + u)) 

e(3,3) = 2i([p 1 p 3 ](-yz + yu + z 2 ) - \pi,p 3 ]u(y + z) 
+ [P2,P3}(y 2 ~yz- zufj/(?,u 2 (y - z + u) 2 ) 

e(4, 1) = 2i([p,p 1 ](-sz -xu + yt) + [p,p 2 ]{-xu + yt) + [p,p 3 ]z(-t + u) 
+ [Pi,P2]s(x + t) + [pi,p 3 ](-xu + yt) 
+ [P2,P3](sz - xu + ytfj/(9tu(x + z + t)(s — t + u)) 

e(4, 2) = -4i(\p, Pl ]s(t -u) + [pi,p 2 ]s 2 + \pi,p 3 ]s(t - u))/{9tu{s - t + u) 2 ) 

e(4, 3) = 2i([p, pi](— xu + yt — zt + tu) — 2[p,p 2 ]xu 

+ \p,P3](sz — xu — zt + zu) + [pi,p 2 ](sy — sz + su + xu) 
+\p2,P3](sy -xu + yt- yu)^/(3u 2 (s -t + u)(y -z + u)) 
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e(4,4) = -2i(jp,p!](st - su - t 2 ) + [p 1 ,p 2 ](s 2 - st - tu) 

-\p 1: p 3 ]u(s + t))/(3u 2 (s - t + u) 2 ) 
e(5, 1) = -4i([p,pi]*(a; + t) + [p,p 2 ]^(^ + - [p,p 3 ]^ 2 )/(9a;t(a; + z + tf) 

e(5,2) = -2i([p,pi](-sz-xu + yt) + [p,p 2 ](-:ra + yt) + [p,Pa]^(-< + «) 

+[pi,p 2 ]s(x + £) + [pi,p 3 ](-xM + yt) 

+ [P2,P3}( s z - xu + ytfj/(9xt(x + z + t)(s — t + u)) 
e(5, 3) = 2i(\jp, pi] (—xz — 2zt + tu) + [p, p 2 ] (s- 2 — xz — xu — zi) 

+ [p,p 3 }(sz -xu- yt) + \p 1 ,p 2 \(-sz + xu- yt) 

+ [Pi, ??,](- + xu + yz + zu) 

+ [P2,P3}(-xy + zu)^j / (3xu(x + z + t)(y — z + u)) 

e(5, 4) = — 2i(jp, pi](— sz — xu + yt) + [p,p 2 ](— xt — xu + yt — zt) 
+ [p,P3](-xt - 2zt + zu) + [p 1 ,p 2 ](sx - tu) 

+ \Pl,P3](-St - XU 

+yt - tu) + [p2,P3,]{sz - xu + ytf)/(3xu(x + z + t)(s - t + u)) 

e (M) = -2i(\p,pi]s(t-u) + [p,p 2 ]s(x-y) + [p u p 3 }s(t - u) 

+ [P2,P3}s(x - y)^/(9xt(s - t + u)(s -x + y)) 

e(6, 3) = — 2i(jp, pi](— sz + xu — yt) + [p,p 2 ](sy — xy + xu — yt) 
+\p,P3](-xz + yu) + \p 1 ,p 2 \(-2sy - su + xy) 
+\pi,P3](xu + yz-yt- yu) 

+ [P2,P3]{sz + xu- ytfj/(3xu(y - z + u)(s -x + y)) 

e(6,4) = —2i(\p,pi](—su + xt) + [p,p 2 }(sx - sy + sz - xu) 
+ [p,p 3 ](sz - xu + yt) + [p 1 ,p 2 ](-sz + xu + yt) 
+ [Pi,P3](-sz + st - su + xu) 

+ [P2,P3}(sx - 2sy - yuf)/(3xu(s -t + u)(s -x + y)) 

e(7,3) = 4i( - [p,P3]z 2 + \pi,p 3 ]z(y + u) + [p 2 ,p 3 ]z(y + u))/(9yu(y - z + uf ) 

e(7, 4) = -2i([p,pi](-zu + yt) + [p,p 2 ](-sz - xu + yt) 

+ [p,P3]z(-t + u) + [pi,p 2 ]s(y + u) + [pi,p 3 ](sz - xu + yt) 
+\P2,P3\(-xu + ytj)/(9yu(s - t + u)(y - z + u)) 

e(8,4) = -2i([p,pi]s(* - u) + [p,p 2 ]s(x - y) + [pi,p 3 ]s(* - M ) 

+[P2,P 3 ]s(z - y))/(9yu(s - t + u)(s -x + y)). 
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The remaining matrix elements can be obtained by the interchange of 
Pi Vii t ^ x , u ^ U from the above results. 



e(5,5) <- 


-> e(l,l) 


e(7,5) <- 


— e(3,l) 


e(8,5) <- 


— e(4,l) 


e(8,8) <- 


— e(4,4) 


e(8,l) <- 


— e(5,4) 


e(8,3) <- 


— e(7,4) 



e(6,5) <- 


— e(2,l) 


e(7,6) <- 


— e(3,2) 


e(8,6) - 


— e(4,2) 


e(6,l) <- 


— e(5,2) 


e(7,2) <- 


— »• e(6,3) 



e(6,6) <- 


— e(2,2) 


e(7,7) - 


— e(3,3) 


e(8,7) <- 


— e(4,3) 


e(7,l) <- 


— »• e(5,3) 


e(8,2) <- 


- e(6,4) 
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Figure captions 

Fig. 1 Diagrams for 7*g — > Gq (a) and 7*G — > (b). 

Fig. 2 Diagrams for 7*g — >■ gGG (a), 7*G — > ggG (b) and 7*g — > qqq (c). 

Fig. 3 Dependence of the various 0(a s ) (a-c) and 0(a 2 ) (e-f) processes on the 
total cms energy y/s. 

a, d) parton densities given by set a: eq. (j33|). 

b, e) parton densities given by set b: eq. (|3 



c,f) parton densities given by set c: eq. fl35|) . 

The Q 2 range is restricted to 5 GeV 2 < Q 2 < s and the cut-off values in 
(0) are y c = 0.04 and M = 2 GeV. 

Fig. 4 x-dependence of the 0(a s ) 3-jet processes at s/s = 20 GeV for each sub- 
process and the sum (solid line) of all contributing subprocesses. 

a) parton densities given by set a: eq. (|33l). 

b) parton densities given by set b: eq. (134]) . 

c) parton densities given by set c: eq. fl3"5|) . 



The Q 2 range is restricted to 5 GeV 2 < Q 2 < 7 GeV 2 and the cut-off values 
in (0) are y c = 0.04 and M = 2 GeV. 

Fig. 5 same as Fig. 4 for 9 GeV 2 < Q 2 < 11 GeV 2 . 

Fig. 6 ^-dependence of the 0(a s ) 3-jet processes at yfs = 20 GeV for each sub- 
process and the sum (solid line) of all contributing subprocesses. 

a) parton densities given by set a: eq. ([33]). 

b) parton densities given by set b: eq. (|34|). 

c) parton densities given by set c: eq. (|35|). 



The Q 2 range is restricted to 5 GeV 2 < Q 2 < 7 GeV 2 and the cut-off values 
in (0) are y c = 0.04 and M = 2 GeV. 

Fig. 7 same as Fig. 6 for 9 GeV 2 < Q 2 < 11 GeV 2 . 

Fig. 8 same as fig. 4 for the 0{a 2 s ) 4-jet processes. 

Fig. 9 same as fig. 6 for the 0(a 2 ) 4-jet processes. 
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